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Abstract. Much work has gone into the construction of quasicontinuum energies that reduce 
the coupling error along the interface between atomistic and continuum regions. The largest 
consistency errors are typically pointwise 0(|) errors, and in some cases this has been reduced 
Gh to pointwise Oil) errors. In this paper we show that one cannot create a coupling method using 

a finite-range coupling interface that has o(l)-consistency in the interface, and we use this to 
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give an upper bound on the order of convergence in discrete w 1,p -norms in ID 

OS 
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1. Introduction 



Atomistic to continuum coupling schemes attempt to provide the high accuracy of an atom- 
istic simulation while also providing the computational savings of a coarse-grained continuum 
approximation. Such simulations are used in the study of localized defects that interact with a 
long-range elastic field, such as crack tips, voids, and interstitials. The quasicontinuum method 
directly couples an atomistic region, where each atom interacts with every other atom within 
a cut-off radius, to a continuum region, whose strain energy density is computed by assuming 
Cauchy-Born kinematics (locally uniform deformation gradients). Many variants of the scheme 
Q\ \ exist that differ on how to couple the regions together. 

One can classify quasicontinuum methods based on whether they create a total energy or 
directly couple forces in such a way that there does not exist a total energy. The original formu- 
lation of the quasicontinuum method |21|. [15] creates a total energy by defining an energy for each 
degree of freedom (atom or finite element) and using a weighted sum. Unfortunately, the model 
contains spurious forces at the atomistic to continuum interface called ghost forces [T7J. These 
forces are a pointwise O(^) consistency error in the linearized scheme and reduce the order of 
convergence with respect to the interatomic spacing parameter e [lj. Subsequent methods such 
as the quasinonlocal quasicontinuum method [20] or the geometrical reconstruction scheme [7] 



removed some ghost forces, but were restricted on the number of neighbors or to specific inter- 
face geometry in 2D and 3D. Recently, the ECC and ACC quasicontinuum variants have been 
proposed in ID and 2D , and these schemes are free of ghost forces for the case of pair-potential 
interactions |18j . A related method, derived as an extension of QNL, is formulated in ID in 

An alternative approach to consistency is the use of the force-based quasicontinuum method [21 
El E], which assigns forces to any degree of freedom, atomistic or continuum, as though the 
entire model were of the same type. Since the continuum model and atomistic model are (^In- 
consistent and there are no special force laws in the atomistic-to-continuum interface, this creates 
a globally point-wise consistent scheme for any geometry. However, this method is more difficult 
to analyze since it does not derive from an energy [H [6] . Solutions to the force-based equations 
can be approximated by employing the ghost-force correction scheme [191 [3] . 

Here we ask the question of whether it is possible to create a quasicontinuum energy which 
is globally o(l)-pointwise consistent. We will find in Section [3] that even in the case of a linear, 
1-D model it is impossible to make an o(l) consistent blending of the atomistic and continuum 
region for any finite-sized interface region. While such a finiteness assumption is common, a 
quasicontinuum force coupling with growing interfaces is treated in |13j . In Section [H we will 
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show how this leads to a bound on the maximal order of convergence in discrete ti; 1,p -norms in 
ID. Convergence analysis for related quasicontinuum schemes include \12\ [8| U] I10j. 

2. Models 

We denote the positions of a periodic atomistic chain by 

y£R N = {R l : y i+N = y e + F for all t G Z}, 

where F denotes the macroscopic deformation gradient for the periodic length and N G £>o 
denotes the number of atoms in the periodic cell. We scale the reference length of the the 
interatomic spacing by e = i in order to have a well-defined continuum limit. We denote the 
displacements from the uniformly deformed position by m = y% — Fei for all i e Z. We define 
the backward difference quotient 

Ui — Ui— r ^7>0 

D r Ui = r G £ , 

re 

where we will write Dui to denote D\Ui. We will write the vector of such differences as Du. We 
likewise define the centered second-difference quotient 

Dl Ui = u i+r ~ 2ui + Ui _ r r g z>0 ^ 

2.1. Atomistic and continuum models. We consider the atomistic energy composed of pair- 
wise interactions with the nearest R neighbors, given by 

R N 
r=l i=l 

where (ft denotes the interatomic interaction potential and R denotes the finite-range cut-off. 
The energy is scaled by the lattice spacing e so that there is a finite energy in the limit N — > oo. 
The continuum approximation uses the same pair potential and approximates the total energy 
by 

R N 
r=l i=l 

That is, interactions of the form 4>{ ~ ^' r ^j are replaced by (^^r ^' - ^'" 1 J . The continuum energy 
is an accurate approximation when yi is smooth. In practical applications, the cost of computing 
the continuum energy is reduced by computing the energy as a sum over the nodes in a piecewise 
linear mesh rather than a sum over all atoms. 
Expanding the pairwise interactions, we have 

<f>{rF + D rUi ) rj <p(rF) + <ft'(rF)D r Ui + (rF){D rUi ) 2 + ■■■ for r G Z >0 . 

For the atomistic model, the (scaled) forces on the atoms are 

{L%u)i = -I^M = -f2r 2 cft"(rF)D 2 r u u (3) 

r=l 

and for the continuum model, the (scaled) forces on the nodes are 



{L c £ u) % = = -J2r 2 cft'\rF)D\ l . (4) 
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2.2. Quasicontinuum coupling. A quasicontinuum energy couples the atomistic and contin- 
uum energies, by partitioning the chain into an atomistic region A and a continuum region C 
that satisfy 

iUC=(0,l] and A D C = 0. (5) 

We assume that A is a finite union of intervals. An atom % is in the atomistic region if si G A. 
Away from the boundary of these regions, the energy is computed using the energy contributions 
£f or £ f. Near the interface between the regions, there can exist special interfacial energies that 
we assume are finite range and which do not change in form with e. Note that the interfacial 
interactions may be longer-ranged than the atomistic ones. 

The above assumptions imply that away from the atomistic to continuum interface and for 
sufficiently large N, the atoms in the continuum region only feel continuum forces (Jl|) and the 
atoms in the atomistic region only feel atomistic forces ([3j). We define A £ and C £ as the interior 
portions of the atomistic and continuum region that only feel interactions from within their 
specific region. We define the interface region I £ = [0, 1] \ {A £ U C e ). In the following, we will 
focus on a single interface between the atomistic and continuum regions. 

3. Consistency at the Interface 
We say that a family of operators L £ is o(l)-consistent with the atomistic operator L a if 

lim II L £ u — Vtu I = (6) 

for every u £ C 2 er (Q), where for every N, u G C 2 er (Q) defines a vector u € ~R. N by U{ = u(i/N). 
Equivalently, we have the local consistency requirement that the equations are consistent if 

(L £ u — L a £ u) = for the three vectors Uj = l,j,j 2 . (7) 

We note that the continuum operator L c e is consistent with the atomistic operator. In the 
following, we show that one cannot satisfy these consistency equations for a quasicontinuum 
energy even in the simplified case of second-neighbor interactions in a 1-D chain with harmonic 
pair potentials. 

We now restrict ourselves to the second-neighbor {R = 2) case for the atomistic and contin- 
uum energies. We write the quasicontinuum operator in terms of its first and second- neighbor 
contributions 

Lf=^{4>"{F)Lf + </>"(2F)Lf), (8) 

where the operators L qc and L| c are independent of (j) and e. This form is possible when we 
assume that the corresponding energy is a function of the "strains" D r u only. We then consider 
possible choices for quasicontinuum energies. Since the first neighbor term for both atomistic 
and continuum operator are identical, it is standard to take the same operator for the first 
neighbor term of a quasicontinuum energy. This gives {Lfu)j = -e 2 D 2 Uj . The second neighbor 
term is a symmetric operator satisfying 

'— Uj+2 + 2uj — Uj_2 = — 4e 2 £)|uj ei € A £ , 

(Lfu)i = I 9 9 (9) 

2 I -Auj+i + 8uj - Auj-i = -Ae 2 D 2 Ui eieC £ , 

< 

where we have not yet specified the interaction law in the interface I £ . We now consider one 
boundary between atomistic and continuum regions. Let m denote the number of atoms in 
the interval of the atomistic to continuum interface I £ surrounding the boundary between the 
regions. This is represented in the operator L q 2 c by an m x m set of coefficients, (L^ij for 
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1 < i,j < m. We have labeled the interface starting at j = 1 and will assume without loss of 
generality that (L| c u)j is continuum for j < 1 and atomistic for j > m. 

Proposition 3.1. There is no symmetric linear quasicontinuum operator Lf c satisfying (|7|). 

Proof. We suppose toward contradiction that L qc satisfies ©• Then, for all i, we have that 

^(4 c -L»-)i = and ^(L^-L«)j 2 = 0, (10) 
3 3 
where j runs over all indices for our periodic chain, though due to the limited range of interactions 
outside the interface, we may restrict j = — 1, . . . , m + 2. Now, we will sum the equations above 
for i = l,...,m, pre-multiplying by i 2 and i, respectively, in order to cancel out the interface 
variables in the interior. We note that by symmetry, L? c = L^- for j < 1 and L q( j = L?- for 
j > m. We have 

m m+2 m m+2 

j=l j=—l i=l j=— 1 

mm m 

= EE(^-^)(^-^ 2 ) + E E (^-w^-v 2 ) 

i=l j=l j=l JG{-1,0} 

m 

+ E E (^-Wj-w 3 ) 

i=l jG{m+l,m+2} 

We cancel the first summation by the symmetry of both L| c and L a , and the third summation 
cancels identically. The second summation has a single non-zero term i = l,j = — 1, so that we 
have 

m 

° = E E (L^-L^j-if) 
i=1 ie{-T0} 
= 1(-1 - 1) = -2. 

We have arrived at a contradiction to our assumption that L qc was consistent (jlOp . □ 

4. Global convergence error 

The result above shows that any linear quasicontinuum scheme must have at least 0(1) 
consistency error in the £°° norm. We can use the above bound on interfacial consistency error 
to give an upper bound on the order of convergence for any quasicontinuum energy. We recall 
the definitions of the discrete if- norms on W N , 

l/p 

\u\\fP := I e y \uj\ p \ , 1 < p < oo, 




ii Uo := max u,- 

E l<j<N 

Provided that Lf c is shift-invariant, that is, it satisfies ([7]) for Vj = 1, we can write the operator 
in terms of differences, 

Lfu = LfDu, 

where we note that the coefficients of eLf are bounded independently from e. We have the 
following lower bound for the error. 



THERE IS NO POINTWISE CONSISTENT QUASICONTINUUM ENERGY 



5 



Corollary 4.1. There exists u £ Cp er (Q) such that for L% c u qc = L®u the error e = u — u qc 
satisfies 

\\De\\ £ P > Ce l+l ' p 

for 1 < p < oo. 

Proof. Since the coefficients of eLf c are bounded, we have the inequality ||L' 3C u||^ CX) < 2 \\Dv\ 

<2 

per 



e 

Using the fact that the operator is not o(l)-consistent, we choose u £ C^ er such that 



We then write 



lim\\(L a £ -L?)u\\ ir >c. 



c < \\(L a £ - Lf)u\\ lr = \\Lfe\\ iT < — \\De\ 



Which implies H-De^oo > Ce. The estimate for other norms follows from the equivalence esti- 
mates for £P norms, \\De\\ e v > e 1 /* \\De\\ l? > Ce 1+1 / p . □ 

This result implies that the rates of convergence obtained in ID for the quasinonlocal qua- 
sicontinuum method in [H |16] are asymptotically optimal among all possible quasicontinuum 
energies with finite-range interface interactions. 

5. Conclusion 

There is much past and ongoing work to create consistent atomistic-to-continuum coupling 
schemes, and a major obstacle are the 0(^)-ghost forces common in such models. Quasicon- 
tinuum energies that are 0( Inconsistent have been constructed in ID and 2D, at least for the 
case of pair-potential interactions. In this work we have shown the impossibility of improving 
the asymptotic truncation error for a linear quasicontinuum energy in ID. Since a higher di- 
mensional problems can behave like a ID configuration due to the geometry (such as simple 
uniaxial strain applied parallel to a planar atomistic-to-continuum interface), this implies that 
it is impossible to make a general o(l)-consistent coupling energy in higher dimensions as well. 
The effect on the rate of convergence shows that in ID the quasinonlocal quasicontinuum energy 
and its generalizations are asymptotically optimal among quasicontinuum energies with respect 
to -u^'P-norm convergence. However, for other quantities of interest, such as the critical strain 
for dislocation movement [SJ, such a method may not be asymptotically optimal. 

Acknowledgements 

The author wishes to thank Mitchell Luskin for many helpful comments on the manuscript. 
This work is supported in part by the NSF Mathematical Sciences Postdoc Research Fellowship 
and ANR PARMAT (ANR-06-CIS6-006). 

References 

[1] A. Abdulle, P. Lin, and A. V. Shapeev. Homogenization-based analysis of quasicontinuum method for complex 
crystals. larXiv: 1006.03781 

[2] W. Curtin and R. Miller. Atomistic/continuum coupling in computational materials science. Modell. Simul. 

Mater. Set. Eng., 11(3):R33-R68, 2003. 
[3] M. Dobson and M. Luskin. Analysis of a force-based quasicontinuum method. M2AN Math. Model. Numer. 

Anal, 42(1):113-139, 2008. 

[4] M. Dobson and M. Luskin. An optimal order error analysis of the one-dimensional quasicontinuum approxi- 
mation. SI AM J. Numer. Anal, 47(4):2455-2475, 2009. 



THERE IS NO POINTWISE CONSISTENT QUASICONTINUUM ENERGY 



6 



[5] M. Dobson, M. Luskin, and C. Ortner. Stability, instability, and error of the force-based quasicontinuum 

approximation. Archive for Rational Mechanics and Analysis, 197(l):179-202, 2010. 
[6] M. Dobson, C. Ortner, and A. V. Shapeev. The spectrum of the force-based quasicontinuum operator for a 

homogeneous periodic chain, submitted. 
[7] W. E, J. Lu, and J. Yang. Uniform accuracy of the quasicontinuum method. Phys. Rev. B, 74:214115, 2006. 
[8] M. Gunzburger and Y. Zhang. Quadrature-rule type approximations to the quasicontinuum method for short 

and long-range interatomic interactions. Comput. Methods Appl. Mech. Engrg., 199:648-659, 2010. 
[9] B. V. Koten, X. H. Li, M. Luskin, and C. Ortner. A computational and theoretical investigation of the 

accuracy of quasicontinuum methods. In I. Graham, T. Hou, O. Lakkis, and R. Scheichl, editors, Numerical 

Analysis of Multiscale Problems, to appear. 
[10] B. V. Koten and M. Luskin. Analysis of energy-based blended quasicontinuum approximations. SIAM J. 

Numer. Anal, to appear. 

[11] X. H. Li and M. Luskin. A generalized quasi-nonlocal atomistic-to-continuum coupling method with finite 
range interaction. IMA Journal of Numerical Analysis, to appear. arXiv.org: 1007.2336. 

[12] P. Lin. Convergence analysis of a quasi-continuum approximation for a two-dimensional material. SIAM J. 
Numer. Anal, 45(l):313-332, 2007. 

[13] J. Lu and P. Ming. Convergence of a force-based hybrid method for atomistic and continuum models in three 
dimension. arXiv:1102.2523 

[14] C. Makridakis, C. Ortner, and E. Siili. Stress-based atomistic/continuum coupling: a new variant of the 

quasicontinuum approximation, preprint. 
[15] R. Miller and E. Tadmor. The quasicontinuum method: Overview, applications and current directions. J. 

Comput. Aided Mater. Des., 9(3):203-239, 2002. 
[16] P. Ming and J. Z. Yang. Analysis of a one-dimensional nonlocal quasi-continuum method. Multiscale Model. 

Simul, 7(4): 1838-1875, 2009. 
[17] D. Rodney and R. Phillips. Structure and strength of dislocation junctions: An atomic level analysis. Phys. 

Rev. Lett, 82(8):1704-1707, Feb 1999. 
[18] A. V. Shapeev. Consistent energy-based atomistic/continuum coupling for two-body potentials in Id and 2d. 

Multiscale Model. Simul., to appear. arXiv:1010.0512 
[19] V. Shenoy, R. Miller, E. Tadmor, D. Rodney, R. Phillips, and M. Ortiz. An adaptive finite element approach 

to atomic-scale mechanics — the quasicontinuum method. J. Mech. Phys. Solids, 47(3):611-642, March 1999. 
[20] T. Shimokawa, J. Mortensen, J. Schiotz, and K. Jacobsen. Matching conditions in the quasicontinuum 

method: Removal of the error introduced at the interface between the coarse-grained and fully atomistic 

regions. Phys. Rev. B, 69(21) :214104, 2004. 
[21] E. Tadmor, M. Ortiz, and R. Phillips. Quasicontinuum analysis of defects in solids. Phil. Mag. A, 73(6):1529- 

1563, 1996. 

CERMICS - ENPC, 6 et 8 avenue Blaise Pascal, Cit Descartes - Champs sur Marne, 77455 Marne 
la Valle Cedex 2 (FRANCE) 



